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CALCULATIONS OF TEE SUPERSONIC WAVE DRAG OF NONLIFTING WINGS 
.WITH ARBITRARY SWEEPBACK AND ASPECT RATIO 

WINGS SWEPT BEHIND THE MACH LINES 
By Sidney M. Harmon and Margaret D. Swanson 


SUMMARY 


On the basis of a recently developed theory for finite swept - 
hack wings at supersonic speeds, calculations of tile supersonic 
wave drag at zero lift were made for a series of wings having thin 
symmetrical biconvex sections with untapered plan forms and various 
angles of sweepback and aspect ratios. The results are presented 
in a unified form so that a single chart permits the direct deter- 
mination of, the wave drag for this family of airfoils for an exten- 
sive range of aspect ratio and sweepback angle for stream Mach 
numbers up to a value corresponding to that at which the Mach line 
coincides with the wing leading edge. 

The calculations showed that in general the wave-drag coeffi- 
cient decreased with increasing sweepback. At Mach numbers for 
which the Mach lines are appreciably ahead of the wing leading ^ 
edge, the wave-drag coefficient decreased to an important extent * 
with increases in aspect ratio or slenderness ratio. At Mach 
numbers for which the Mach lines approach the wing leading edge 
(Mach numbers approaching a value equal to the secant of the angle 
of sweepback), the wave-drag coefficient decreased with reductions 
in aspect ratio or slenderness ratio. In order to check the results 
obtained by the theory, a comparison was made with the results of 
tests at the Langley Memorial Aeronautical Laboratory of swept - 
back wing attached to a froely falling body. The variation of 
the drag with Mach number and aspect ratio as given by the theory 
appeared to be in reasonable agreement with experiments for the 
range of Mach number tested. 


INTRODUCTION 


Recent developments in airfoil theory for supersonic speeds 
(references 1 and 2) indicate pronounced effects of sweepback and 
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aspect ratio on the drag. In reference 1, a theory was developed 
for calculating the pressure distribution at supersonic speeds and 
zero lift for swept -back wings of arbitrary linear taper and aspect 
ratio having thin symmetrical airfoil sections of sharp leading 
edges . 

In the present paper, the method of reference 1 is applied to 
calculate the supersonic -wave drag for a series of wings having 
thin symmetrical biconvex sections at zero lift with untapered 
plan forms and various angles of sweepbaclc and aspect ratios. The 
term "biconvex profile" as used herein refers to an airfoil section 
composed of two parabolic arcs. In each case, the wing is con- 
sidered to be cut off in a direction parallel to the direction of 
flight. In the calculations the Mach number is varied from 1.0 
to a value corresponding to that at which the Mach line coincides 
with the wing leading edge. The results of the calculations are 
presented in a unified form which permits the direct determination 
from a single chart of the wave drag for this family of wings for 
an extensive range of sweepback angle .and aspect ratio for Mach 
numbers from 1.0 to the value corresponding to that at which the 
Mach line coincides with the wing leading edge, or equal to the 
secant of the angle of sweepback. . Although the calculations have 
been made for the biconvex profile, the data may be applied to 
indicate corresponding results for profiles approximately similar 
to the biconvex. 


In order to illustrate the possible applicability of the 
theory to a typical swept-back wing, the calculated dra g of a wing 
of sweepback at zero lift is compared with the results of 
tests made at supersonic speeds at the Langley Memorial Aeronautical 
Laboratory on swept-back wings attached to a freely falling body. 


SYMBOLS 


x, y, z 

dz/dx 

c 

t/c 

A 


coordinates of mutually perpendicular system of axis 
in wing 

slope of airfoil surface 

chord of airfoil section, measured in flight direction 
thickness ratio of section, measured in flight direction 
angle of sweep, degrees 


m = cot A 
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h wing semispan measured along y-axis, semichords except 

in appendix A 

K parameter indicating spamd.se position equal to y/m, 

semi chords 


> 


A 

l/t 

P 

p/<i 

V 

u 


u 

V 


I 


aspect ratio 

slenderness ratio, ratio of wing semispan measured along 
leading edge to maximum thickness of center section 

disturbance pressure 

pressure coefficient, ratio of disturbance pressure to 
dynamic pressure in free stream 

velocity in flight direction 

x-component of disturbance velocity, positive in flight 
direction 

u caused by source line with reversal in sign of m 

z- component of disturbance velocity 

disturbance-velocity potential 

source factor required to maintain a given wedge angle 


M Mach number 



y coordinate measured along y-axis which is shifted to tip 

meeticr, 'semicbprds 

y coordinate measured along y-axis which is shifted to 

opposite tip section, semi chords 

d wave drag at section 



wave-drag coefficient at section without tip effect 


c d 

^d 


wave-drag coefficient at section including tip effect 

increment in section wave-drag coefficient caused by wing 
tips # 4 



k 
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Ac^ 


-II 


nD 

c 

C D 

AC. 


D 


AC D-, 


AC35 


II 


^Dp 

£> n 


increment in section wave-drag coefficient caused ’ey wing 
tip located on same half of wing as section 

increment in wave-drag coefficient at section on one wing • 
panel caused by tip of opposite wing panel 

wave -drag coefficient for wing without tip effect 

wave-drag coefficient for mng including tip effect 

increment in wave-drag coefficient caused by tips, com- 
plete wing 

increment in rave-drag coefficient on one wing panel caused 
by adjacent wing tip. complete wing 

increment in rave-drag coefficient on one wing panel 
caused by tip of opposite wing panel, complete wing 

drag coefficient obtained as sum of coefficients for rave 
drag and friction drag 

auxiliary variables which replace x and y , respec- 
tively, used to indicate origin of source line 


Primed values of A, y , y Q . y^, h, m, dz/dx, and z indi- 
cate transformation involving multiplication ey factor P 


Subscripts: 

1. 2 wings related according to transformation which makes yP 
and mp eq.ua 1 respectively for two wings. 

Subscript notations for u and u indicate tne origin of source 
line in terms of coordinates x and y, respectively 

ANALYSIS 


Basic data .- The present analysis is based on thin-airfoil 
theory for small pressure disturbances and a constant velocity of 
sound throughout the fluid. The axes used are the mutually per- 
pendicular x, y, z system in which the x-axis is taken in the 
direction of flight positive to the rear, the y-axis is along the 


ft 
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span positive to the right, and the z-axLs is positive up-wards. 
Figure 1 shows the symbols used to designate the wing-plan-form 
parameters. The present analysis is made for untapered brings of 
biconvex profile at zero lift and is limited to a Mach n um ber 
range from 1.0 to the value corresponding to that at which the Mach 
line coincides with the wing leading edge, that is, to a value 
equal to the secant of the angle of sweepback. 

Theory . - If p is the disturbance pressure computed for one 
surface of the airfoil section, the wave drag for the section is 


d = 


2 


1 c 
l*0 



( 1 ) 


and the section wave-drag coefficient 


c d 



c j q. dx 

HO 


(£) 


where dz/dx is the slope of the surface of the airfoil at the 
point x. For the symmetrical biconvex profiles (composed of 
parabolic arcs) 



where the thickness ratio t/c may be considered the sole shape 
parameter. From thin-airfoil theory, 


P _ _ £ ^ _ gU 
q. ~ " V dx " V 
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where u is the disturbance- velocity component in the x-direction, 
taken positive in the flight direction. For a given svept-back 

wing with A = cot \i, the drag coefficient in equation (2) at 
the spanwise station y may now be written as 


c d (y) 





m fC 


y/m 




(3) 


The desired integrand u in equation (3) may be determined 
by the procedure given in reference 1. On the basis of the 
linearized theory, reference 1 derives a solution representing an 
oblique (svept-back) source line making the angle of sveepback A 
with the y-axis. This solution utilized for the pressure field 
or for the disturbance velocity is the real part of 


Uq, 0 = I cosh 


x - mp 2 y 
P.|y - mx j 


( 4 ) 


where the subscript notation indicates that the 
at the origin of coordinates (x = 0, y = 0) * 
shown in' reference 1 to satisfy the boundary co 

/dzY 

oblique wedge making the ha If -angle rn 

coordinate system of reference 1 (y' = yP, s* 


source line starts 
Equation (4) is 
ndition for a thin 

the transformed 
= zp) where 


(~\ = — - - AiY ,. 1. in order to obtain an equal wedge angle 

\dx j V Y m' 

in the physical coordinate system, the following relations between 
the transformed coordinates of reference 1 and the physical coor- 
dinates are used 


m 1 = mp 

*. 2 - 

hz ' 

V = Sz 
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vhere w' and v refer, respectively, to the vertical velocities 
in the transformed and physical coordinate systems. Thus 


w 


50 dz ' 

dz ' dz. 


w'p 


The half -angle of the wedge is then determined as 

dz w _ jt \|l~- 
dx ~ V " V m 

or the required source factor in order to maintain the desired 
wedge angle is 


I 


V m dz 


(5) 


By superposition of solutions of the wedge type, swept -"back 
wings of desired profile shape and plan form can "be "built up 
(reference 1) . In order to satisfy the "boundary conditions over 
the surface of a "biconvex wing, semi-infinite source lines of 
equal strength are placed along the leading and trailing edges 
"beginning at the center section, in conjunction with a constant 
distribution of sink lines along the chord. At the tip, where the 
wing is cut off, reversed 3emi-infinite source and sink lines are 
distributed so as to cancel exactly the effect of those originating 
at the center, section in the entire region of space outboard of the 
tip. In the present analysis, the tip is assumed to be cut off in 
the direction of flight;- The term "tip effect” refers to the 
effect of this wing cut-off . The form of the integrand u for 
equation ( 3 ) is given in appendix A. 

In calculating the wave drag over the wing, the disturbances 
due to the elementary sources, and sinks are evidently limited to 
the regions enclosed by their Mach cones. Figure 2 shows the 
typical Mach lines originating from the source lines at the leading 
and trailing edges of the center and tip sections over a wing of 
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45° sweepback. Figure 2(a) shows the Mach lines for the infinitely 
long wing, and figure 2(b) includes the Mach lines starting from 
the tip section. In each case the disturbance over the wing Caused 
by the leading- and trailing-edge source lines is limited to the 
region of the wing behind the corresponding Mach line. The regions 
affected by each of the Mach lines are indicated in figure 2(b) as 
regions I to IF. Region I represents the part of the wing affected 
by the sink line starting from the leading edge of the adjacent 
tip; region II represents the wing area affected by the sink line 
starting from the opposite tip. Region III is influenced by the 
source line starting at the leading edge of the center section and 
includes the entire wing; region IF is influenced by the source 
line starting at the trailing edge of the center section. 


The resultant velocity u at a point on the wing is made of 
the component velocities caused by each of these source and sink 
line 3 where the influence of each component is restricted. to the 
region behind its Mach line. The drag coefficient C & c;) there- 
fore obtained by evaluating in equation (3) the integrand Uq >0 
over the entire section (region III), the integrand u c ^q over 

part of the section included in region IF, and the integrands for 
the u-components caused by the sink distributions along the pro. 
file (fig.. 2(b) and appendix A, equation (A2)). The. drag coeffi- 
cients Ac,* and Ac d __ are obtained similarly by integrating 
^*1 XX 

along the section in the regions I and II, respectively, in addi- 
tion to the integrations for the u-components caused. by the source 
distributions along the profile (appendix A). The limits of inte- 
grations for x along the chord and for y along the span, which 
represent the boundaries for the regions of influence for tho 
individual u-functions required for a biconvex profile, arc given 
in tho table in appendix A. 


Formulas for section wave-dr ag coefficients .- The formulas 
for tho section drag coefficients obtained by integration of the 
u-functions and by use of equation (3) are presented. in appendix 
These formulas give expressions for the drag coefficient without 
the tip effect c d and also the expressions for the increments 

in c j caused by the tip effect Ac d . 


Wave -drag coefficients for complete wing .- In the present 

investigation the section 'dra'g coefficients tho 

equations in appendix B wore integrated graphically to obtain t 
results for tho wing- drag coefficients. Subsequent y, > 
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analytical expressions for the integrations were obtained. These 
formulas for the wing-drag coefficients are presented in appendix B. 

Drag coefficient of swept-back wing 'at Mach number of 1.0 .- 
The solution of the equations for c^ given in appendix B shows 

that, for a symmetrical untapered finite swept-back wing at Mach 
number of 1.0 and zero lift, positive and negative infinite values 
for c^ are obtained at various sections of the wing. The inte- 
gration over the wing of the limiting values for these infinite 
terms, however, gives zero. Although some sections of the wing 
have infinitely positive or negative drag, the total drag coeffi- 
cient over the ling results in a finite value. The prediction of . 
infinite values of drag at certain sections of the wing clearly 
violates at these sections the assumption of small disturbances 
from which the linearized theory is derived. The calculated values 
for the total drag coefficient at Mach number 1.0 are therefore 
questionable. The formulas for the total drag at Mach number 1.0 
are presented in appendix B. 

Conversion of drag s olution to series of related wings .- An 
examination of equations (B.1) , (B3) and (B5) indicates that the 
drag solution obtained for one value of m and M may be applied 
directly to obtain the drag for a -whole series of -wings in which 
each wing is at a certain. appropriate Mach number. (Equation (Bl) 
is formed by adding expression (Bib) to the right-hand side of 
equation (Bla).) For example, equation (Bl) may be expressed in 
the following form: 



where F (*&> mp\ refers to the variable terms and wrhere m£ = m', 

W / 

and - = K. If the subscript 1 refers to a wing at Mach number 
m 

corresponding to P 1 and the subscript 2 to any other wing at the 
Mach n umb er corresponding to then the drag coefficients for 

the two wings may be obtained from equation (6) as followrs. 



8 /t \ 2 m I p l 
* \c ) 1 Pi 




(7a) 
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(7b) 


Equation (7) shows that if = y 2 P 2 


and 



(t/c) 1 - Pg 


= c d 

(t/c) 




(t/c )} 2 cot 
(t/c) cot Ag 



where Cj and Cj refer to the spamd.se positions y-, and 

u w 2 

respectively. In a similar manner, it can be shown that if 
p 2 . 

two wings are related according to y^P^ = y 2 p2 ail & m iP.l “ m 2^2> 

the section drag coefficients obtained for wings 1 and 2 from equa- 
tions (B3) and (B$) are in the same ratio as that expressed in 
equation (8), Equation (8), therefore, may be generalized to apply 
to the total drag coefficient at the section or, 


c d]_ * c do 


(t/c)! 2 cot A]_ 
. (t/c) 2 2 cot A 2 


( 9 ) 


where c^ and C£ 0 refer to the spanwise positions y^ and 

yiPi 

, respectively . 

P2 


The wing-drag coefficients for wings 1 and 2 are given, 
respectively, by 


'•Pi bh 


n h. 


uo 


; d x d 71 


( 10 ) 
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Cd 2 h 2 


A he 


uo 


: d 2 d ?2 


( 11 ) 


By substituting for the integrand c 4 in equation (10) the rela- 

dy 2 p 2 


tionship expressed in equation (9) and by substituting 
for dy x , equation (10) may be written as: 


Pi 




(t/c)^ cot Aq 


0h c 


h 2 (t/c) 2 ' cot. Ag Jo 


c d 2 


or 



(t/c)^ cot A x 
(t/c) 2 2 cot A 2 


( 12 ) 


Equation (12) permits a rapid determination oF^bhe'drag^coeffi- 
cient for wings of arbitrary sveepback, aspect ratio, and thickness 
ratio (within limitations of airfoil theory) from data obtained for 
one svept-back wing for the appropriate range of aspect ratio and 
Mach number. For this purppse, use of a wing, of 45° sveepback as 
the reference wing is most convenient. If the subscript 2 is used 
to refer to the parameters for the wing of 1+5° sveepback and the 
subscript 1 is dropped, equation (12) becomes 


c D 2 (V°)^ cot A 

(t/o) 2 a . 


( 13 ) 


where Cp and Cp refer to wings whose aspect ratios and Mach 
£ 

numbers are related by the following equations: 
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A 2 = A tan 

A 

(14) 

P 2 = P cot 

A ■ 

(15) 


The foregoing analysis shows that the results obtained for a 
•wing of 45° sweepback and a. given aspect ratio A 2 can he trans- 
formed hy means of equation ( 13 ) to all wings for which the aspect- 
ratio parameter A tan A = A 2 and the Mach number parameter 

P cot A = p£- 

The grouping of the parameters as indicated in the foregoing 
analysis permits the use of a single generalized chart for pre- 
senting the drag results.. This chart is discussed in the section 
entitled 'Results and Discussion. " 

Prandtl-Glauert rule modi fi ed for sup er soni c flow s . - In con- 
sidering the linearized problem of supersonic flow past a wing, it 
is often convenient to refer the supersonic results for a given 
wing to a transformed wing at a reference Mach number of M- = 2. 

If this transformation is used a rule resembling the Prandtl-Glauert 
rule (reference 3) for the subsonic case, where M = 0 is 
the reference Mach number, may be. obtained. .This rule may be stated 
as follows: 

The streamline field of the supersonic flow for a given body 
at a stream Mach number M may be calculated by multiplying the 
given y- and z- dimensions, including those for the Mach linos, by 

the factor \/M“ - 1 and then by calculating .the flow about the 

resulting transformed body at the Mach number' V 2. The pressure p 
and velocity increments u for the given body at the Mach number M 
can then be obtained by multiplying the calculated pressure p and 
velocity increments u at corresponding points of the transformed 

body by the factor — 

M 2 - 1 

It is interesting to note that the derivation of formulas (13) 
to ( 15 ) as given in this paper corresponds to utilizing the solution 
for a transformed wing for the whole family of wings related to 
this transformation and then applying the aforementioned modified 
Prandtl-Glauert rale. 


! 
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RESULTS AND DISCUSSION 


Variation of section drag coefficient along span .- Figure 2, 
which was introduced previously to illustrate the system of Mach 
lines, al30 shows the variation of section drag coefficient c d 

along the Bpan. The data are presented for a wing of 45° sweepback 
and thickness ratio of 0.10. Figure 2(a) gives the results for a 
wing of infinite aspect ratio and figure 2(h) gives the results for 
a wing of finite span. 

• In figure 2(a), the data are shown for Mach numbers of 1.100, 
1'3^3, and 1.4l4. The lowest Mach number 1.100 represents a case 
in which the wing leading edge diverges rapidly from the Mach line 
(upper part of fig. 2(a)). In this case, the section drag coeffi- 
cient c d has a maximum value of 0.0542 at the center section, 
then drops sharply to zero at a distance of 1.13 chords from the 
center line . Beyond this point, the wing shows a negative drag, 
which approaches asymptotically the subsonic value of zero at an 
infinite distance from the wing center. . This type of ■wave-drag 
distribution is similar to that indicated in figure 11 of refer- 
ence 1 for a wing of 60 ° sweepback at a Mach number of 1.4. 

For the higher Mach number 1.343, the spamri.se variation of c d 

is markedly flatter. Unlike the preceding case, the drag coeffi- 
cient does . not have. its maximum value at the center section but 
at first, increases in the outboard direction, then reaches a peak 
and falls ter zero -at- a distance- from -the- center _of_ 6.6 chords 
(not shown in fig. 2(a)). ' " 

At highest Mach number 1.4l4, the Mach line becomes coin- 
cident with the wing loading edge. In this caso, the wing gives 
a very high drag and the section drag coefficient increases in the 
outboard direction, approaching infinity at an infinite distance 
from the wing center. 


Figure 2(b) illust rates the condition at which the aspect ratio 

is less than 1 A/m 2 - 1 . The calculated case shown is for an aspect 
ratio of 1.86 and Mach number of 1.10. In this case, the two wing 
tips cause increments in section drag coefficient on each half of 


. o uu uaun lid jl 

the wing, namely, Ac dl and Ac d . The tip effect Ac dT at 


given distance from the tip is independent of the aspect ratio, 
tip effect Ac d ^_, however, is a function of the aspect ratio. 


The 
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Effect of wing tip on wing-drag coefficient *- Figure 3 shows 
the typical variation with aspect ratio of the increment in Cp 
due to the tip ACp. The data are shorn for a wing of 45° sweep- 

hack, and. a thickness ratio of 0.10 for Mach numbers of 1.100, 1*343} 
and 1.414. The present analysis for the untapered wings indicated 
that if the aspect ratio is equal to or greater than 2m/ (mp + l), 
the integrated value of Ac^j over the wing is zero. On this basis, 

if the aspect ratio of the wing is greater than i/\|m^ - 1, the 
total increment in drag contributed by tip is zero. As the aspect 
ratio is reduced, the tip effect ACp^_ which occurs when the 

aspect ratio i3 smaller than lj\fu^ - 1, however, leads to an 
increase in Cp . The tip effect ACp then reaches a peak value at 

a certain aspect ratio, but as this aspect ratio is further decreased, 
ACp drops sharply. In this case, at aspect ratios of approxi- 
mately 0.5, ACp becomes zero and as3iunes large negative values, 
with further reductions in aspect ratio. For. applications to very 
small aspect ratios, however, the theory may require' modifications.' 
The data, in figure 3 show that as the Mach number is increased, the 
aspect ratio corresponding to zero value of ACp becomes smaller. 

The tip effects shown in figure 3 for the wing of 45° sweep- 
back are similar for other wings of different swoepback at appro- 
priate aspect ratios. The conversion formulas indicated in the 
section entitled "Analysis" indicate that the aspect-ratio effects 
for wings of different sweepback correspond qualitatively for equal- 
values of the aspect -ratio. parameter A tan A* Tho Mach, numbers 
for each of the wings differ, however. An aspect ratio of 0.8. 

(fig. 3) for the wring of 45° sweepback at a Mach number of 1.10 
for example corresponds to an aspect ratio of 0.8 cot A at a Mach 

number equal to \jl~ + j(l.l0)^ - lj tan^ A for any other wing of 
sweepback angle A • 

Variation of wing-drag coefficient with Mach n umber , sweepback, 
and slenderness ratio .- Figure 4 shows the variation of Cp with M‘ 

for different sweepback angles with constant slenderness ratios. 

The slenderness ratio represents the ratio of the wing semispan 
measured along the leading edge to the maximum, thickness of the. 
center section. The data, are presented for sweepback angles of 30°, 
45°, 52.5°, and 60° with slenderness ratios of 25 and 50. The 
wings in figure 4 for tho different slenderness ratios and sweep- 
back angles are assumed to have the same wing area and the same" 
profile normal to the wing loading edge. The slenderness ratios 
aro based on a thickness ratio of 0.10 measured in a direction 
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normal to the wing leading edge. The thicloieBS ratio . t/c measured 
in' the flight direction, therefore,, varies with sweepback as cos A 
or is equal to 0.1 cos A. • The aspect ratio is reduced with sweep- 

back by the factor cos 2 A. The aspect ratio is related to the 
slenderness ratio. by the formulas 

; A = 0.2 ^1^cos 2 A 

The plan .forms for the different wings, are shown in figure 4. 

The results in figure 4 show that, in general, the drag coef- 
ficient decreases with increasing sweepback. At Mach numbers for 
which the Mach lines are appreciably ahead of the wing leading .edge, 
increasing the slenderness ratio or aspect ratio gives important 
reductions in calculated wave-drag coefficient. At Mach numbers for 
which the Mach lines' approach the wing leading edge (M— + sec A), 
however, short wide wings give appreciable reductions in wave-drag 
coefficient. The figure also indicates that the highest wave-drag 
coefficients for the normal range of aspect ratio occur at a Mach 
number equal to sec A. 


Effect of aspect ratio on wing-drag coefficient .- Figure 5 
indicates the effect of aspect ratio on the wave-drag coefficient 
for the wing. The data in figure 5 show the wave -drag- coefficient 
Cp tan A 

parameter — — ~ plotted against the aspect -ratio parameter 

A tan A These re sults a re shown for~various values of~thcrMach — 

number parameter \/m - 1 cot ,A, which correspond to a range of 
Mach number from 1.0 to the secant of the angle of sweepback, 
or 1 < M < sec A. . . .... , 


Figure 5 shows that for a given value of the Mach n umb er 
parameter, the maximum wave-drag coefficient occurs at a definite 

aspect ratio. For example, if A = 45° and ^M 2 - 1 cot A = 0.310 
(that is, M = 1.05) the maximum value of Cj» occurs at an aspect 
ratio of 0.05, if the aspect ratio is decreased to values sma ller 
than 0.05, C D drops very sharply. Similarly, as the aspect ratio 

is increased from 0 .85, also decreases. Thus, in general for 

ij b« . h i 

the Mach dumber, parameter corresponding to \Jk 2 — 1 cot A = 0.310 
the maximum value of occurs at an aspect ratio equal to 
O.85 cot A 

of c urves of figure 5 to wings of arbitrary sweep - 
back an d aspect ratio . - The scale labels and mirvAH n-p fs j irc ^ 
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apply to the series of wings that may be derived from a basic 
10 '■percent -thick, 45° swept -back wing. The labels express the 
transformation equations (13) to (l^). If A is set equal to k?° 
and t/c is set equal to 0.10 in these labels to correspond to 
the basic wing, the ordinates become simply Op, the abscissa A 

(aspect ratio), and the curve parameter \/m^ - 1. The results in 
figure 5 may be applied to all swept-back wings covering a range 
of aspect ratio from.. 0 to 10 cot A corresponding to a range of 
Mach number from 1.0 to sec A* The data apply specifically to 
untapered Tangs with biconvex profiles at zero lift with the wing 
tips cut off in the direction of flight. The results, however, may 
be applied to indicate approximate results for profiles similar to 
the biconvex. . 

The following example is given in order to illustrate the use 
of figure 5. For a given wing 


A = 70° 

A. - 3 

- = 0.08 
c 

M = 2.20 • ‘ 


In order to find Cp*. 


A tan A = 8.2k 


\/m 2 - 


1 cot A 


0.715 


From figure 5, for A tan A = 8.2k and \/M 2 - 1 cot A - O.715 


C-n tan A 

— = 0.0123 

100 (t/c y- 
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therefore 


c D = 0.00286 


Comparison of theory with experimental results . - Figure . 6 
shows a comparison of the theoretical results for wave-drag coeffi- 
cient with data obtained from unpublished tests made at the Langley 
Laboratory for swept -back wings attached to a freely falling body. 

The comparison is shown for a wing of 1 - 5 ° sweepback for two aspect 
ratios at zero lift. The aspect ratios for the experimental wings 
based on over-all span were 3.6 and 5 * 4 . Because of the small 
thickness of the wing section relative to the' thickness of the ' 
fuselage (O.72 in. to 10.75 in.), the calculations assumed. that 
the fuselage acted as a perfect reflection plane. On this basis, 
the aspect ratios of 3*6 and 5*4 corresponded to aspect ratios 
of 2.65 and 4 . 42 . The experimental wing had an NACA 65-OO9 profile 
taken in the direction normal to the wing leading edge or a thick- 
ness ratio of O.O636 in the flight direction; whereas, the cal- 
culated results are based on a biconvex having the s am e ma ximum 
thickness ratio as the experimental profile. The calculated drag 
coefficient Cp was obtained by adding the skin-friction drag 

coefficient of 0.0027 to the calculated wave-drag coefficient. This 
value of O.OO27 for the skin-friction drag coefficient was obtained 
as the minimum drag coefficient from the unpublished experimental 

results and this value appears somewhat, lo w. The comparison is 

shown for a range of Mach number from 1.0 to 1.257 The plan form 

of the wings are shorn in figure 6 . 

The comparison in figure 6 indicates that the calculated 
values of the wave-drag coefficients are of the order of magnitude 
of the experimental results • The agreement appears to be closer 
for the wing of higher aspect ratio than for the wing of smaller aspect 
ratio* It should be noted that in addition to the approximations 
inherent in the linearized theory, the calculations do not take 
account of factors such as fuselage interference, boundary-layer 
effects, and the exact profile. 


CONCLUDING EEMARKS 


A theoretical investigation has been made of the supersonic 
wave. drag of untapered swept-back wings at zero lift. The wing 
section 'was biconvex and the wing tip wa s considered to be cut°off 
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in the direction of flight. The investigation -was limited to a. 
range of stream Mach number from 1.0 to a value corresponding to 
that at which the Mach line coincided with the wing leading edge. 

For this range of Mach number, the following conclusions have been 
draw.: 1 

1. In general, the calculated wave-drag coefficient decreased 
with increasing sveepback. 

2. At Mach numbers for which the Mach lines are appreciably 
ahead of the wing leading edge, increasing the slenderness .ratio 
or aspect ratio gave important reduction in the calculated wave- 
drag coefficient. 

3* At Mach numbers for which the Mach lines approach the wing 
leading edge (Mach numbers approaching a value equal to the secant 
of the angle of sveepback), decreasing the slenderness ratio or 
aspect ratio reduced the calculated wave-drag coefficient. 

4. The highest calculated wave-drag coefficients for the 
normal range of aspect ratio occurred at a Mach number equal to 
the secant of the angle of sweepback. 

5 • The maximum wave-drag coefficient occurred at a definite 
aspect ratio which is determined by the Mach number and angle of 
sweepback. 

6. For aspect ratios greater than i/\/m 2 - 1, where M is 
the Mach number, the increment in wave -drag coefficient .for the 
wi.ng contributed by the tip was zero. 

7 • The variation of the drag with Mach number obtained for 
one sweepback angle for appropriate aspect ratios may be presented 
in a unified form so that the drag for the complete range of sw r eep- 
back angle, aspect ratio, and Mach number may be directly determined 
from a single chart. 

8. A comparison of the theory for a Mach number range from 1.0 
to 1.29 with results obtained from tests of swept -back wings attached 
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to a freely falling "body indicated that the calculated values were 
of the order of magnitude of the experimental results. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya. 
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APPENDIX A 


F0EMU14S POE THE INTEGRAND u IN EQUATION (3) FOE FINITE 


UNTAPERED SWEPT -BACK WINGS OF BICONVEX PROFILE 


AT EERO LUT 


m ss cot A < — 
~ P 


In order to satisfy the boundary condition for a finite swept- 
back wing of biconvex profile, the integrand u in equation (3), 
may be expressed in terms of components caused by elementary source 
.lines as follows: 



where the subscript notation indicates the origin of the source 
line. The bars over u refer to the soiree lines caused by the 
opposite wing panel; that is, u indicates a source line with a 
reversal in the sign of m. 

In equation (Ai), the u-functions are given by the real parts 
of the following expressions: 


u 


y* to* 


(x. y) = I cosh 


-!*■£“ m|3 2 (yr - tj) 


P y ■ r a(x - g) 


u (x, y) a I cosh — 

S ' 11 B|y - 


+ mp 2 (y - t>) 
tj + m(x - g) | 


P 


hi p 
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where rj represents the origin of the elementary source lilies. 
For the biconvex profile. 


= ot 

dx C c 


and from equation (5), the factor 


1 = 2— — 

c Jt 


m 


\/l - m 2 p 2 


The symbol jjj in equation (Al) refers to an integration opera 

tion which represents the influence of the uniform distribution of 
sink lines along the chord of the biconvex profile. This symbol 
is d efined by the following expressions: 


„(x, y) = 

O 7 ! 


■>o 


Cx-py 


dl . -1 
— cash 


x - g - mp 2 (y - ti) 

Ply ■ n " m(x - g) ! 




= (y - n) 


n — 

dl j \/l - m 2 p 
dx I m 


cosh” 1 


* - s 

P )y - tj| 



U 


m(x - g) 




cosh 


- mp 

-1 P(y - _n) 


ll - 


»(■ 


L “ | ) 

• T> 


(A2) 


y - P 
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where, for a biconvex profile 


dl _ V m d^z 
dx * to a ' 

= ht/o) l — a — 

C ” \jx - m 2 V 2 


Equation (A£) may be expressed as a function of that io 

y - a 


V** y> - (7 - ,) r(^) 


Then 


Au ( x, v) 


*(y - ft) f 


iL 1 


; (y - ti) J 


The limits of integrations with regard to x for the section 
drag coefficients and with regard to y for the total wing-drag 
coefficients are discussed. The u-dbmponents caused by each of the 
elementary source linos are aero at all points outside of the 
respective Mach cones. The functions for the u- integrand in equa- 
tion (3) are therefore evaluated along the section for values of x 
beginning at the forward boundary of the Mach cone. This integra- 
tion gives the section-drag-coefficient components. Similarly in 
order to obtain the wing-drag coefficient, the section drag coeffi- 
cient components obtained from the respective u-functions are 
evaluated along the wing span for values of y contained within 
tne Mach cone. The following table refers to one side of the wing 
[x and y positive) and shows the limits of integration for x 
and y for the required u-functions: 


WlM £ I «IH 
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Limits of integration 


u- components 

X 

y 

j Lower 

1 limit 

Upper 

limit 

Lower 

limit 

Upper 

limit 

U 0,0 ^0,0 
l" 0,0 ' 5^0,0 

y/m 

£ + c 

. m 

0 

h 

U C ,0 u c, 0 - 

1 .1-. 

D U c,0 p u c,0 

yP + c 

z + c 
m 

0 

h 

i 

u h/m,h 

1 

D U h/m,h 

1 

g(mg + 1 ) - y (3 

c ' 

. me 

la •• 

mp -i- 1 

(if A > 2 m ^ 

i h 


XU 

\ mp 4 * 1 / 

0 (if A < : 2m \ 

V “mp + 1/ 


u . 

h/m, -h 

1- 

-h 

— — i 


2+ c 
m 

_L 

h(mp -i- 1) - me 
1 - mp 

(if A > - £m ^ 

\ mp + 1 / 

0 fit A < 2m \ 

\ “ mP + 1/ j 

1 

L 

h 



J 
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FORMULAS FOE WAVE -DB AG COEFFICIENTS FOE FINITE UNTAPEEED 
SV/EPT-BACK WINGS OF SYMMETEICAL BICONVEX PROFILE 
AT ZERO LIFT 


m = cot A < -• | 

- P J 


Section Drag Coefficients 


In the following analysis the quantities y and K are 
employed nondimensionally in terms of the semichord. The equations 
for the drag coefficients in all cases refer to the real parts of 
the indicated expressions. 


Section drag coefficients without tip effect .- The section drag 
coefficient for the given wing at a apamri.se position y and Mach - 
number M without the tip effect was 'found to be as follows: 


The term 


2m 


represents a convenient integration limit 


1 - m0 

which indicates the intersection of the Mach line from the center- 

section trailing edge with the wing leading edge. For y = Km < ^ 

1 - mp. 

C d <y> . 8 (if M f oosh -l KjlS . 2 ooah -i m 

» « Vc/ y \ Km' »'/ 

L k(i - a' 2 ) + 2 


+ — <.2 cosh 


3M 


2m 1 


fl - m' 2 1. v 

2^2K^ - 3K - l) cosh 

4K(2K S - 3 ) cosh" 1 B^|sc(jL - m’ 2 ) 

cia 1 1 


:;w3 . w . i i “I 01 + 2 

2m 1 (K + 1) 


\/(l - m ,£ ) (K + 2) 2 - (Km') 2 j 




(Bla) 
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For y = Km > 
tion (31a) . 


2m 

1 - TB&’ 


the following expression is added to eq.ua- 


8/t ^ Jk 3 co 3 h- 1K 


-(-) 
n v « c / 


3\jl - m’ 2 i- 


Km' 
2 cosh"^ 


-1 k(i - in’ 2 ) 


2m 1 


2 ( 2IC 3 - 3K + l) cosh" 1 


2m'(K - 1) 



where m' = mp. 


(Bib) 


For the special case m ■== — 


m 2 - 1 

mth the wing leading edge, and the expression for the drag coeffi- 
cient obtained as a limiting case (m } — > 1) for all values of y 
becomes: 


the Mach line coincides 


CD 


V 

OJ 

f>\ 

y' 3 coalf 1 y' + 2 - £ 

y'( 3y' + h) (y ' - 1 ) - 2 

\oJ 

y' 3 

~ ” > 

\/y * + 1 


(B2) 


where 


y' = yP 


At the center section, where y or K = 0, 



equation (Bla) becomes: 
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3« 



m 



cosh 


1 I_ 

m3 


At the center section, for y = 0 and m = -, equation (B2) 
becomes: 



Increment in section drag coefficient caus ed by wing tips .- 
The increment in c^ caused by the tip3 depends on various factors, 

such as the sweep angle, aspect ratio, and Mach number. The following 
types occur in an untapered wing: 

I> i Y , the aspect rati0 of the wing is ec.ua 1 to or greater 
than l/\/M 2 - 1 each tip affects solely its own half of the wing. 

In this case the effect of the tip is limited to the region of the 
wing outboard and rearward of the front Mach line originating from 
this tip. (See. fig. 2(b).) The region of the wing affected°is 

between values of y from h - — -- — to h. 

m3 + 1 

The increment in section drag coefficient at a Mach number M 
and spanwise position y caused by the tip was found to be as 
follows: 


Acgfy) = - ei ^ 

-*• 7C \ C / 


ya 


12m' 


(, 


,2 2 + m 1 


m ' 


- 12m ' ) cosh 


iV + 2m ' / 

i f Pm » 

- 37a) 


,2 

*r ! cm 

m 'y * V 

a 


" y a 


cosh 


•1 (l - m' 2 )y • + 


2m’ 




B1 


,2 


(B3) 
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vhere the subscript a indicates that the x-axis is shifted to the 
tip section, and y a ' = y Q P and m' = m3 . In the plan form of the 

•wing 


T = b + 


In equation (B3) values for y a may he taken from - 


cm 


m3 + 1 

to 0. When the Mach line is coincident with the leading edge of the 
airfoilj that is, m = — — t the expression for Ac^ becomes: 


a t \ ^ /t \ 2 

Ac d j(y) = - (-) m 

ic '-c / 


\/m 2 - 1 

V(V 2 - ■*) 


cosh 


-1 V + 2 


IV 


\fV 


&y a ,£ - 2y a ' + 4 ) 


(B4) 


II* If the aspect ratio of- the wing is less than — — — , the 

\/m 2 •• 1 

tip on the opposite wing contributes an increment in c, in addi- 

d 

t ion “to'-that _ dfscus sed~under _ type— Ii — The _ ±ncr ement~ in — at~a 

section caused by the opposite tip was obtained in the following form: 


Ac, 


dn v 


<x> *!(;/ “ wC 77 * ’ ' 10h ’ ' ' ■ 2Ch ' -«’>] 2 - (»’jv •)* 

^ 1 


.2 


(l4 + m'~) - 3h 'y, ' + 2h ,2 -m' 2 


12 


cosh 


-1 y b ' " £ (b' -m') 

~V . 


3m^l - m' 2 U 


2 y h |3 - Gy h '^h ' + 3y, ' (2h - m' 2 ) - 2h ’3 + 3m ,2 h ' - m 


cosh 


, ± yb't + *' d ) - 2 (b ' - m') 

£m'^y b ' - h' + m'^ 


(B5) 
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where the subscript b indicates that tho x-axis is shifted to the 
opposite tip section, and where y b ' a y^P, m' = nip, gri d h' = hp. 

In the plan form of the tang 


y = y b - h. 


The limits for y^ to he used in equation (B5) depend on the 
value of aspect ratio A relative to the parameter . . . Thus 


mp + 1 

(a) If the aspect ratio of the wing is greater than 


2m 


mp + 1 

the front Mach line from the opposite tip intersects the trailing 
edge at a value of y b = ~ ( - h - so that values for y in 

U 1 - -mfl J h 


equation (B5) may he taken from — — to 2h at the tip. 


(See fig. 2(h) . ) 


mp 


(b) If the aspect ratio A is equal to or less than 


2m 


mp + 1 

the front Mach line from the opposite tip intersects the center 
section and values for y^ in equation (B5) may be taken from h 

to 2h. In this case,, the increment in Ac^ discussed under 
type I is obtained at spanwise positions of y g from -h to 0. 


When m = 


{F- 


Z> equation (B5) becomes: 


" 7 (1 ) n \ h ' h ( 7y *> ' ' 10h ’ ' ‘ 2 ) yV ' ~ 2 ( h : ' - 1 ) 2 . - y b 

v* L 

- (%b “ 3h \ ' + 2h ,2 - 1 j (cosh 


,2 



2 ! 2 Jb ’ 6y b ,dh ' + 3y b ' (sh ,2 - l) 

+ — 
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The total Increment in drag coefficient at a section caused by 
the tips is given by 


Ac, = A ca + Acj 
d a.j T 


and the total drag coefficient at the section is 


c a = 


c d 


CD 


Ac, 


Wing-Drag Coefficients 


Wing-drag coefficient without tip effects.- If the aspect ratio 


is equal to or less than 
the tip effect is 


2m 


1 - OP 


■> the wing-drag coefficient without 




A' 2 ( , -l A' + 2m' 

1 3A ' cosh — 


12m' 3 


A 'm' 


6a' cosh 


~ - UA' 2 (l - m’ 2 ) + 4m ' (A ' + m’) 
m' v 


+ 2A 


\|l - m 


• 2 + i. 


3m 


°v 1 - 


nr 


a*' 3 oosh -i Aji .'1 + an' 

,2 


2m' 


(&a' 3 + 3A'm' S - A' 3 ) cosh' 1 *1^ t m '^ + £m ' 

2m'(A' + m*) 


(2A’ 3 - 6A’m' 2 ) 


cosh 


-1 1 + m 


,2 


where A' = AP and m' = mp 
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If the aspect ratio is greater than 
coefficient -without tip effect is, 


2m. 


1 - m| 3 


the wing-drag 


<v!GJ 


8/tf m f A" 


,12m 


,3 


, ^ cr _,.-l A' + 2m' -1 A' - 2m '\ 


3A ' I cosh 

\ A 'm ' 


+ cosn 


A 'm * J 


- 6a' cosh -1 ~ - \Ja ' ^ (l - m'") + 4m ’(A' + m') 
m' - 


- \iA ,l "(l - m'^) + 4m'(m' - A') + 2A'\|l - m 


-X A'(l - m' j + 2m' 
cosh ’ — 


3m ,3 \jl - m' ; 


r 

t *2 

<2m' 3 


o> , 

1C ! 

i 

J 


p 

2m 1 


- cosh 


-1 A'(l - m' 2 } - 2m' 


am'" 


(3A'm' 2 - 2m' 3 - A' 3 ) cash' 1 

2m '(A 1 - m') 


,2 0 ,3 , ,3V . -l.A'\l + n'v + 2 m' 

'm ■+ 2m' - k' J j cosh — — — 

2m ' (A ' + m') 


( 0 . .3 , ,2^ ,-11 + m' 2 

+ l,2A - oA m'./ cosh — 

2m’ 

Increment in wing-drag coefficien t cau s ed by w ing tips.- If 
the aspect ratio of the wing is equal to or greater than l/|3, the 
contribution of the wing tips to the wing-drag coefficient is zero. 
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If the aspect ratio of the wing is equal to or greater than 


2m 


the total increment in Cp caused by the local tip, or Z£p^, is 

zero,* and the increment ACp is obtained solely from the effect of 
the tip located on the opposite half of the wing, or ADp^. For 
this case, integration of equation (B5) over the wing yields: 

f X 


m3 + 1 




2m , )(A* + m') 2 cosh" 1 - + 1 

A* + m' 


+ (A 1 + 2m') (A 1 - m') 2 cosh" 1 A SJi 

|A ' - m’( 


A 1 

3m'3 


6m 


,2 


A' 2 (2 + m' 2 )] cosh" 1 - — Jl - A’ 2 i 
-» A* 3m'" I 


If the aspect ratio of the wing is less than 


2m 


the 


mp + 1* 

increment ACp is affected by both wing tips. For this case, 
integration of equations (B3) and (B5) yields: 




m 


v 3m’3^l - m' 2 


(A* - 2m')(A* + m')^ cosh 


-1 m'A 1 +"1 
A’ + m' 


+ (A' + 2m*)(A' - m’) 2 cosh' 


1 1 a,' 3 coeh” 1 a«:-A 

|A’ - &’} 


2m 


,2 


- (a»'3 - 3m ,2 A' + A-3) cosh' 1 ^1— ± 


{.2m'(m' - A’) } 


+ ^n 

i 

A' 3 


6m' 2 - A '^(2 + m ,2 /j cosh" 1 l- + ^4m'(m' - A’) +A‘ 2 (l -m ,2 )J 


4m 


,3 


coeh-* ^ - A ’ - 


m’A ' 


3m 1 


1 - A 


,2 
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Total wing-drag coefficient .- The total wing .drag is obtained 
as the sura. 


On = Cj) + ACp 

CO 


where the components jand are calculated from- the fore- 

going equations for the wing-drag coefficient appropriate to the 
aspect ratio of the wing. 

Wing-drag coefficient for special case m = When the Mach 
£ 


line coincides .with the wing leading edge /m = -V the wing- drag 

. - \ . fi l - v . : 

coefficient obtained for all aspect ratios as a limiting case (m 1 — } 1) 

is equal to the real part of the following expression: 



Wing-Drag Coefficients at Mach Number of 1.0 

The drag coefficient for the wing at:; M = 1.0 - may be expressed 
in terms of the following formulas which were obtained by integrating 
over the wing the limiting values for c^ in equations (Bl), (B3), 

and (B5) as the factor p = \j M 2 - 1 approaches zeros 
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than 


°D = 


where 


(l) If the aspect ratio of the wing is equal to or greater 
2 m 


3* 


ft 2 


8 (c) m|(-K 3 + 6 K + k) log e (IC + 2) - 3K 3 log e 2K 

\!/0 


(4K 3 - 6K - 2) log (K + 1) + - 2) 


dK 


ph/m 


r. 

(-IC 3 + 6K - 4) log (K - 2 ) - (k 3 - 6K - 4) log (K + 2) 


(4E 3 - 6 k) log (k 2 - l) + 2 log — — - 6K 3 log It] 

e e K + 1 6 J 


d IK 


,*h-2m 


(y ~ h) 3 + 3(h - y) 


2m- 


m 


log 0 2(y - h + 2m)J 


(y + h ) 3 , 3 y(y - h ) . 3 h + y ) r , . 

2m 3 m 3 m J " 3 

% 3 6 y A 


.m- 


3 m 


" | log Q 2m(y + m) j 


f Z._l2i(7y - 3h - 2m)(y - h + 2m) > dy' 


4m 3 


J 


(B6) 


K = — 
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34 


MCA EM No. L 6 K 29 


In equation (b 6) , the first two integrals represent the drag 
coefficient for the wing without the tip effect; whereas, the last 
integral represents the effect of the tips For this case, 

in which the aspect ratio is equal to or greater than . 2m, the tip 
effect is zero; hence the integral for the section increments 

in c^ is not given in equation (B6) . Equation (B6) has been 

solved for a sweepback angle of 45° and the results for this sweep 
angle may be converted to other wings of arbitrary sweepback by the 
formulas (13) to (15) • For the wing of 45° sweepback, m = 1 and 
and A > 2, equation (B6) yields the following result: 


^ 1 


(-A 3 + 12A + 1 &) log Q (A + 2) 


+ (“A 3 + 12A - 16) log Q (A - 2) + 2 a(a 2 - 12 ) log Q A - 4A 


(B7) 


(2) If the aspect ratio of the wing is smaller than 2m, the 
upper limit of the first integral in equation (Bb) is reduced from 2 
to h/m, the second integral vanishes, and the lower limit for the 
third integral is reduced from h - 2m to zero*' For this case, 
however, in which the aspect ratio is less than 2m, ZsC^ is not 

zero, and the following integral must be added to those in equa- 
tion (B6) to obtain C^: 


3a W h 


10 


-h 


ya 0 3y a \ 

log 

/2m 3 . e 


!(y a + 2m) 




+ 2 log 


m 


y a + *a 


y a (im - 3y a )(y a + 2m) 


4m- 


hy. 


(38) 


where 


= 


- h 


(/ 
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For the ving of 45° sueepback, m = 1, and A < 2, equations (B 6) 
and (b 8) yield the same result for Cp as that obtained for values 

of A greater than 2, as expressed by equation (B?) . In this case, 
the real part of log e (A -2) is used. 
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